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A polymorphism in the histone deacetylase 1 gene 
is associated with the response to corticosteroids 
in asthmatics 

Min-Hye Kim 1 ' 2 , Sae-Hoon Kim 1 ' 2 ' 3 , Yook-Keun Kim 1 ' 2 '*, Soo-Jong Hong 4 , Kyung-Up Min 1 ' 2 , 
Sang-Heon Cho 1 ' 2 , and Heung-Woo Park 1 ' 2 

Background/Aims: Recent investigations suggest that histone deacetylase 1 
(HDACi) and HDAC2 may be target molecules to predict therapeutic responses to 
corticosteroids. We evaluated the effects of variation in HDACi and HDAC2 on the 
response to corticosteroids in asthmatics. 

Methods: Two single nucleotide polymorphisms (SNPs) were selected after rese- 
quencing HDACi and HDAC2. For the first analysis, we evaluated the association 
between those SNPs and asthma severity in 477 asthmatics. For the second anal- 
ysis, we evaluated the effects of these SNPs on lung function improvements in 
response to corticosteroid treatment in 35 independent adult asthmatics and 70 
childhood asthmatics. 

Results: We found that one SNP in HDACi (^1741981) was significantly related to 
asthma severity in a recessive model (corrected p = 0.036). Adult asthmatics who 
were homozygous for the minor allele of ^1741981 showed significantly lower % 
forced expiratory volume in 1 second (%FEV 1 ) increases in response to systemic 
corticosteroids treatment compared with the heterozygotes or those homozygous 
for the major allele (12.7% ± 7.2% vs. 37.4% ± 33-7%, p = 0.018). Similarly, childhood 
asthmatics who were homozygous for the minor allele of ^1741981 showed sig- 
nificantly lower %FEV 1 increases in response to inhaled corticosteroid treatment 
compared with the heterozygotes or those homozygous for the major allele (14.1% 
± 5.9% vs. 19.4% ± 8.9%, p = 0.035). 

Conclusions: The present study demonstrated that ^1741981 in HDACi was sig- 
nificantly associated with the response to corticosteroid treatment in asthmatics. 
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INTRODUCTION 

Asthma, a chronic airway inflammatory disease, is an 
important source of morbidity and mortality world- 
wide [1]. Current guidelines recommend corticosteroid 
treatment for the management of asthma [2-4]. Endog- 
enous corticosteroid levels and therapeutic responses 



to exogenous corticosteroids are known to be influ- 
enced by genetics, with heritability ranged from 0.40 
to 0.56 [5-7]. This suggests a pharmacogenetic basis 
for the intraindividual variability of corticosteroid re- 
sponsiveness in asthmatics [8]. 

Histone deacetylase (HDAC) plays a key role in the 
regulation of inflammatory genes by removing acetyl 
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groups from histories [9,10]. It has been reported that 
conditional deletion of HDACi in T cells resulted in 
enhanced airway inflammation and increased TI12 cy- 
tokine production [11]. Moreover, HDAC is involved in 
the mechanism of action of corticosteroids. For exam- 
ple, recruitment of HDAC2 to activated inflammatory 
genes is an important mechanism of inflammatory 
gene repression by corticosteroids and HDAC2 activity 
is reduced in some diseases in which patients showed 
poor responses to corticosteroids treatment [12]. To- 
gether, HDACi and HDAC2 may be target molecules 
to predict therapeutic responses to corticosteroids. 
However, to date, any genetic effects of HDACi and 
HDAC2 on the response to corticosteroids in asthmat- 
ics have remained unknown. 

In the present study, we evaluated associations be- 
tween asthma severity, which is related to reduced re- 
sponsiveness to corticosteroids [13], and single nucleo- 
tide polymorphisms (SNPs) in HDACi and HDACz in 
adult asthmatics. Then, we further assessed the effects 
of those SNPs on responsiveness to corticosteroids, as 
measured in terms of the change in lung function, in 
independent adult and childhood asthmatics. 

METHODS 

All subjects enrolled in this study provided written 
informed consent. The study protocol was approved by 
the Institutional Review Board of Seoul National Uni- 
versity Hospital (H-1110-130-384). 

Primary analysis 

In total, 477 asthmatics were enrolled at Seoul Nation- 
al University Hospital, Seoul, Korea. The diagnosis 
of asthma was made at least 6 months prior to enroll- 
ment, according to current guidelines [2-4]: episodic 
symptoms including wheezing, coughing, and dys- 
pnea plus a positive bronchodilator response (an in- 
crease in forced expiratory volume in 1 second [FEVJ, 
from baseline that was more than 200 mL and more 
than 15% of the prebronchodilator value), or a provoc- 
ative concentration of methacholine causing a 20% 
reduction in FEV 1 (PC 2D ) <, 16 mg/mL [2-4]. Enrolled 
asthmatics were treated with conventional medica- 
tions based on the The Global Initiative for Asthma 
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guideline according to their asthma control status [2]. 
Asthma severity was determined based on lung func- 
tion and the medication use index needed to obtain 
control, as described previously [14,15]. 

Secondary analysis 

To generalize the results obtained in the primary 
analysis, we enrolled adult (treatment-naive 35 subjects, 
systemic corticosteroid group) and childhood asthmat- 
ics (treatment-naive 70 subjects, inhaled corticosteroid 
group) in a secondary analysis. Adult asthmatics vis- 
ited our clinic first with dyspnea and decreased lung 
function (FEV-l < 80% predicted value) and were diag- 
nosed with asthma based on positive bronchodilator 
responses. They were treated with short-term systemic 
steroids alone for prompt relief of their symptoms (oral 
prednisolone 15 mg, twice per day for 7 days). Child- 
hood asthmatics enrolled at the Asan Medical Center, 
Seoul, Korea, were diagnosed with asthma according 
to current treatment guidelines [2-4]. All were treated 
with inhaled corticosteroids for 8 weeks according to 
their asthma control status. Treatment responses in 
both groups were measured by the degree of increase 
in FEVj: %FEVj increase = (FEVj after treatment - FEVj 
at baselineyFEVj^ at baseline x 100. 

Genotyping 

The assessment of the SNPs in HDACi and HDAC2 
was described in our previous study [16]. Briefly, after 
isolating genomic DNA from the peripheral blood of 
24 healthy Korean subjects using the QIAamp DNA 
blood kit following the manufacturer's protocol (Qia- 
gen, Hilden, Germany), we amplified 2 kb of the s'-up- 
stream region in the promoter and all exons, includ- 
ing the exon-intron boundaries, of HDACi and HDAC2 
by polymerase chain reaction (PCR; reference genome 
sequences; NM_oo4o64.2 [HDACi] and NM_ooi527.3 
[HDAC2]). Amplified PCR products were sequenced us- 
ing a Big Dye Terminator Cycle Sequencing Ready Re- 
action Kit (Applied Biosystems, Foster City, CA, USA) 
in both directions according to standard protocols. 
After sequencing of HDACi and HDAC2, we identified 
five SNPs in HDACi and 14 SNPs in HDAC2 (Sup- 
plementary Tables 1 and 2). Among these, two SNPs 
(^1741981 in HDACi and ^58677352 in HDAC2) were 
selected for scoring after considering minor allele fre- 
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quencies (higher than 5%) and location (in exons and 
5', near the gene). Scoring was conducted with the high 
throughput single base-pair extension method (SNP- 
IT assay) using an SNPstream25K system, which was 
customized to automatically genotype DNA samples in 
384-well plates and to provide a colorimetric readout 
(Orchid Biosciences, Princeton, NJ, USA) as described 
previously [17]. 

Statistical analyses 

An association analysis, based on a case-control de- 
sign, was performed for each SNP using a genetic 
model approach: an allele model (A [major allele] vs. B 
[minor allele]), a dominant model (AA vs. AB + BB), and 
a recessive model (BB vs. AA + AB). Hardy-Weinberg 
equilibrium was tested using 2x2 tests. All statisti- 
cal analyses were performed with the R version 2.15.3 
(http://www.R-project.org/). p values of less than 0.05 
were considered to indicate statistical significance. 



RESULTS 

The characteristics of the asthmatics enrolled in the 
primary and secondary analyses are shown in Tables 1 
and 2, respectively. RS1741981 in HDACi and ^58677352 
in HDAC2 were in Hardy-Weinberg equilibrium. The 
primary analysis revealed that ^1741981 was signifi- 
cantly related to asthma severity in a recessive model 
(p = 0.006) (Table 3). This relationship was still signif- 
icant after the Bonferroni correction (multiplied by 
6 [2 SNPs x 3 models]; corrected p = 0.036). However, 
^58677352 showed no relationship with asthma sever- 
ity (data not shown). As ^1741981 showed a significant 
relation with asthma severity in a recessive model, the 
secondary analysis was done using the same (reces- 
sive) model. Adult asthmatics with the CC genotype 
of rsi74i98i (n = 12) showed significantly lower %FEV 1 
increases in response to systemic corticosteroid treat- 
ment compared with those with the CT or TT geno- 
type (12.7% ± 7.2% vs. 37.4% ± 33.7%, p = 0.018) (Fig. 1A). 



Table 1. Characteristics of asthmatics enrolled in the primary analysis 



Characteristic 


Mild 


Moderate to severe 


No. 


134 


343 


Male sex 


44(32.8) 


122 (35.7) 


Age, yr 


47.3 ± 17.0 


50.2 ± 16.3 


Basal FEVi predicted, % 


95.2 ± 15.8 


78.6 ± 13.5 


PC20-methacholine a , mg/mL 


6.9 ± 1.2 


4.7 ± 1.1 


Atopy 


66 (49.3) 


157 (45-8) 



Values are presented as number (%) or mean ± SD. 
FEV 17 forced expiratory volume in 1 second. 
a Geometric mean. 



Table 2. Characteristics of asthmatics enrolled in the secondary analysis 



Characteristic 


Adult 


Childhood 




No. 


35 


70 




Male sex 


17(48.5) 


34(48.6) 




Age, yr 


56.7 ± 12.1 


10.2 ± 1.4 




Basal FEVi predicted, % 


57.1 ± 12.9 


78.9 ±7.4 




PC20-methacholine a , mg/mL 


NA 


9.7 ± 1.3 




Atopy 


10 (28.6) 


52 (74-3) 





Values are presented as number (%) or mean ± SD. 

FEVj, forced expiratory volume in 1 second; NA, not applicable. 

a Geometric mean. 
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The same tendency was observed in childhood asth- 
matics. Subjects with the CC genotype of 151741981 (n 
= 15) showed significantly lower %FEV 1 increases in re- 
sponse to inhaled corticosteroid treatment compared 
with those with the CT or TT genotype (14.1% ± 5.9% 
vs. 19.4% ± 8.9%, p = 0.035) (Fig- iB). The combined p 
value, calculated from the one-sided p values of the 
secondary populations using Stouffer z-transform test 
[18], was 0.0019. However, ^58677352 showed no rela- 
tionship with lung function improvement in response 
to corticosteroid treatment adults or childhood asth- 



KJIM^ 

matics (data not shown). 
DISCUSSION 

The present study demonstrated that ^1741981 in 
HDACi is significantly associated with the response to 
corticosteroid treatment in asthmatics. To our knowl- 
edge, this is the first report of a genetic association 
between HDACi and response to corticosteroids in 
asthmatics. 



Table 3. Genotype frequency of 13174198 in HDACi and 1558677352 in HDAC2 among asthmatics in the primary analysis 



Genotype 




Frequency 




p value 




Mild 


Moderate to severe 


Dominant 


Recessive 


Allele 


HDACi rsi74io8 


CC 


17 (12.7) 


83(24.5) 


0.203 


0.006 


0.014 


CT 


75 (56.0) 


171 (50.4) 








TT 


42 (31.3) 


85(25.1) 








HDAC2 ^58677352 


AA 


3 (2.2) 


15 (4-4) 


0.887 


0.406 


1.000 


AG 


51(38.1) 


119 (34-7) 








GG 


80 (59.7) 


209 (60.9) 









Values are presented as number (%). RS174198 genotyping failed in four moderate to severe asthmatics. 



p = 0.018 
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Figure 1. Forced expiratory volume in 1 second (FEVj improvement in response to corticosteroid treatment according to the 
genotype of ^174198 in HDACi in the secondary analysis. (A) Adult asthmatics treated with oral prednisolone 15 mg, twice per 
day for 7 days. (B) Childhood asthmatics treated with inhaled corticosteroids according to their asthma control status. %FEV 1 
increase = (FEV, after treatment - FEVj at baseline) / FEVj at baseline x 100. 
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Corticosteroids are among the most effective treat- 
ments for asthma [19], and reduce inflammation via 
glucocorticoid receptor-mediated recruitment of 
HDAC2 [20]. Because reduced HDAC2 expression was 
observed in bronchial biopsies of asthmatic patients 
[21] and HDAC2-mediated deacetylation of the cortico- 
steroid receptor enabled nuclear factor-KB suppression 
in airway epithelial cell line and alveolar macrophages 
[22], the discussion of the corticosteroid response in 
asthma has to date focused on HDAC2 alone. However, 
several investigators have argued that corticosteroid 
action in asthma is explained partially by inhibition 
of T cell activation [23-25]. Recently, it has been report- 
ed that that T cell-specific loss of HDACi resulted in 
enhanced allergic airway inflammation, by modulat- 
ing cytokine production [11]. Moreover, HDACi was lo- 
calized within most airway cells and infiltrating in- 
flammatory cells in the lung in a murine asthma 
model, and the HDAC inhibitor TSA attenuated aller- 
gic airway inflammation in mice by reducing T cell 
infiltration and TI12 cytokine production [26]. Taken 
together, HDACi may play an important role in medi- 
ating the anti-inflammatory action of corticosteroids. 
Accordingly, our data provide pharmacogenetic evi- 
dence that a genetic polymorphism in HDACi can pre- 
dict the response to corticosteroid treatment in asth- 
matics. 

Functional variants strengthen the reliability of re- 
sults of a pharmacogenetic study. RS1741981 is located 
in the 5' near gene region of HDACi. We used a 
web-based tool to evaluate the functional relevance of 
rsi74i98i. The Functional Single Nucleotide Polymor- 
phism (F-SNP; http://compbio.cs.queensu.ca/F-SNP/) 
database integrates information obtained from 16 bio- 
informatics tools and databases about the functional 
effects of SNPs [27]. For 151741981, F-SNP used TF- 
Search (http://www.cbrc.jp/research/db/TFSEARCH. 
html) [28] and Consite (http://asp.ii. uib.no:8o9o/cgi- 
bin/CONSITE/consite) [29] tools. Both tools computa- 
tionally predict changes in transcription factor bind- 
ing according to the allele of ^1741981 (T allele vs. C 
allele) (Fig. 2 and Supplementary Fig. 1). F-SNP indicat- 
ed 0.268 for the FS score of ^1741981. The FS score is a 
functional SNP scoring system provided by F-SNP to 
distinguish features of disease-related SNPs versus 
neutral SNPs [30]. The median FS score for neutral 
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Figure 2. Changes in transcription factor binding according 
to the allele of ^1741981 in HDACi. (A) -1332T. (B) -1332C. 
Predicted results from Consite (http://asp.ii.uib.no:8oc)o/ 
cgi-bin/CONSITE/consite). 



SNPs is 0.1764, whereas for disease-related SNPs, the 
median rises to 0.5 [30]. Thus, ^1741981 may be func- 
tionally important. 

There are several potential limitations to our re- 
sults. First, sufficient statistical power, conferred by a 
large sample size, is an important aspect of genetic 
studies. The small number of patients enrolled in the 
present study limits the scope of our results. Second, 
subsequent replication studies should consider as 
many clinical and demographic confounding factors 
as possible and should be performed in diverse ethnic 
backgrounds. Third, haplotype approaches for HDACi 
would result in more precise associations. For exam- 
ple, one previous study showed the advantage of exam- 
ining both individual SNPs and haplotypes to further 
pin-point a potential causal SNP [31]. Finally, although 
a computational method suggested that ^1741981 may 
be functionally important, further mechanistic stud- 
ies are needed to clarify the precise role of ^1741981 
and HDACi in the response to corticosteroid treat- 
ment. 
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KEY MESSAGE 

1. In the present study, we identified 151741981 in 
HDACi, which showed a significant relationship 
with asthma severity. 

2. We also found that this genetic variation was 
associated significantly with lung function 
improvements in response to systemic corti- 
costeroid treatment in adult asthmatics and in 
response to inhaled corticosteroids treatment in 
childhood asthmatics. 

3. This is the first report of a significant associ- 
ation between a genetic variation in HDACi 
and the response to corticosteroid treatment in 
asthmatics. 
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